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Abstract 



We study the renormalisation group flows between minimal W models by means of a new 



o 
o 



(N 



set of nonlinear integral equations which provide access to the effective central charge of 
both unitary and nonunitary models. We show that the scaling function associated to the 
■ nonunitary models is a nonmonotonic function of the system size. 

1. A recent study of the renormalisation group flows between nonunitary minimal models 
revealed an unexpected behaviour for the groundstate energy Eq(R), in that it was a non- 
monotonic function of the system size R Q. The nonmonotonicity was illustrated using the 
finite-size scaling function c e g (r), which up to the bulk term is proportional to the groundstate 
^ ; energy 

E (R)=E hulk (M,R)-^^- , r = MR , (1) 

where M is the so-called crossover scale (the mass in massive theories). As the system size 
goes to zero c e fj(r) becomes the effective central charge 

O 1 

£!' lim c eff (r) = c-24A . (2) 

We denote the actual central charge by c while Ao is the conformal dimension of the lowest 
primary field of the UV CFT. 

The effective central charge and the central charge of the unitary minimal models coincide, 
and according to Zamolodchikov's c-theorem @ there exists a function c which is monotonic. 
However, apart from the UV and IR points at which c equals the central charge of the relevant 
CFT, it is not clear if there is any connection with Eo(r). Nevertheless the groundstate 
energy of the unitary models is always monotonic. Analogously it had been thought that 
the groundstate energy of the nonunitary models would also be monotonic, but the results 
of [[jj and |3]-|6| provide a number of counter examples. In this letter we study a further set of 
perturbed conformal field theories, demonstrating that c e fj(r) behaves nonmonotonically for 
the majority of nonunitary models. 

We consider the minimal models WQ^ q based on one of the simply laced Lie algebras 
Q = A n -i, D n , Eq,Et, Eg 0. The models are specified by two coprime integers p and q with 
p > h, in terms of which the central charge and the effective central charge are 

c = N ( l _Hh±2)ii^l) ca = N ( 1 _Mt±}l). (3 , 

V pq ) \ pq J 

Here ./V denotes the rank of the algebra and h the dual Coxeter number. The primary fields 
&n ci' are labelled by a pair of weights f2, W which satisfy 

e-n< q ,e-n'<p, (4) 
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where 9 is the highest root of Q || . All models have a primary field <I> a( jj that has weights = 
fijd and $7 = J7 a dj corresponding to the trivial and adjoint representations of Q respectively. 
It has conformal dimensions 

A = A = 1-^, (5) 

and is relevant for all p, q such that q > p. Formally denoting the action of the unperturbed 
CFT Aqft, that of the perturbed model may be written 



A = A C ft + \ JdTx®^ , (6) 

reproducing for A\ the well-known ^13 perturbations of the Virasoro minimal models. 

Depending on the sign of the coupling constant A, the perturbation either leads to a massive 
quantum field theory, or it induces a 'massless flow' into a conformal field theory with smaller 
effective central charge. One of the standard methods of studying the groundstate energy 
of both types of model is the thermodynamic Bethe ansatz. The result is a set of coupled 
nonlinear integral equations (the TBA equations), whose solution provides direct access to 
c e fj(r) at all values of r. 

The unitary WQ P / +1 can alternatively be described as x Q^/Q 1 ^ 1 coset models at 
k = p—h, in terms of which the perturbing operator <£ a dj is usually known as ^id^adj • 
TBA equations describing the evolution of the effective central charge between these coset 
models, and therefore the unitary WQ models, are already known |H-Fll|, and they verify the 
conjectured pattern of flows |8, 12 1: 



wg p / +1 + $ adj -> w&fr 1 * . (?) 

By dropping the restriction q = p + 1 , we may also consider flows originating from the much 
larger class of nonunitary WQ minimal models, which also have a description as a coset but 
at fractional level Q p ^ q p ^ h x Q^/Q p J ; ^ q p ^ h+1 _ Analogous to the known behaviour of the 
<^13 perturbations of the Virasoro models [13-16], it is natural to suppose that the nonunitary 
generalisation of (0) will be 

wg^ + Q^^wg 2 *-™ . (8) 

TBA equations describing massless flows from these nonunitary models are not yet known. 
Instead, motivated by we propose a different type of nonlinear integral equation (NLIE) 
whose solution provides access to the effective central charge of both unitary and nonunitary 
minimal models. The equations can be found in section 2, and are tested in section 3. In 
4 and 5 we extract some exact results and make a comparison with ultraviolet and infrared 
perturbation theory. The connection between the massive models and the Gross-Neveu models 
and other comments can be found in section 6. 

2. Our starting point is a set of nonlinear integral equations that encode the groundstate 
energy of the imaginary-coupled simply laced affine Toda field theories. Our interest in these 
theories lies in the fact that for values of the Toda coupling constant /3 2 = p/(p + 1) the 
theory can be consistently restricted to the massive ^adg perturbation of the unitary WQ^ P+1 
minimal models [17-p2]|. Moreover, the massive c/>i3 perturbation of the nonunitary WA p ' q 
minimal models can be obtained from the sine-Gordon model by tuning the coupling to 
(3 2 = p/q. A similar result may be true for the nonunitary minimal models based on the other 
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M a = Ms'm(ira/h)/ s'm(ir/h) 1 < a < n - 1 

M n _i = M n = M/2sm(n/h) , M a = M sm(ira/h)/ sm(ir/h) l<a<n-2 
Mi = M 2 = M , M 3 = 2Mcos(vr/4) , M 4 = M 5 = 2Mcos(vr/12) 
M 6 = 4Mcos(vr/12) cos(vr/4) 

Mi = M , M 2 = 2Mcos(5vr/18) , M 3 = 2Mcos(vr/9) , M 4 = 2Mcos(vr/18) 

M 5 = 2M 2 cos(vr/18) , M 6 = 2M 3 cos(2vr/9) , M 7 = 2M 3 cos(vr/18) 

Mi = M , M 2 = 2Mcos(vr/5) , M 3 = 2Mcos(vr/30) , M 4 = 2M 2 cos(7vr/30) 

M 5 = 2M 2 cos(2vr/15) , M 6 = 2M 2 cos(vr/30) 

M 7 = 4M 2 cos(vr/5) cos(7vr/30) , M 8 = 4M 2 cos(vr/5) cos(2vr/15) 

Table 1: The mass normalisation. 



simply laced Lie algebras p3| , and at the level of the NLIEs we do find the choice j3 2 = p/q 
yields both unitary and nonunitary perturbed models. 

The NLIEs describing the groundstate energy of the massive imaginary-coupled Toda 



field theories were first obtained in 24, 25 1, and have appeared in a different context in 
The effective central charge is defined in terms of iV functions which satisfy a set of coupled 
equations 



N 

-{■m a re d + m ^ C^a h + 2*^ [ d9' ip ab (8 - 6')5mln(l + e 
b=i b=i ^ c 



0) 



The integration contour C runs just below the real axis while r is built from the lightest mass 
M of the theory and the cylinder size R via r = MR. We have set m a = M a /M, where 
each mass M a is associated to a node of the Dynkin diagram via the labelling of (27j, and 
is such that (Mi, M 2 , . . . , Mn) forms an eigenvector of the Cartan matrix with eigenvalue 
4sin 2 (7r/2/i). Our particular normalisation is given in table |]. The kernel functions 



>ab 



sinh(%(e+l)fc)) 
sinh(f cosh(f jfe) 



c~ b Hk) 



(10) 



are written in terms of the 'deformed' Cartan matrix C a b(k), which is equal to 2 if a = b, and 
— 1/ cosh(irk/h) if nodes a and b of the relevant Dynkin diagram are connected. Note that at 
k = it reduces to the standard Cartan matrix C a b- The exact effective central charge may 
be determined using 



c e s(r) 



— ^m a d6 sinh 9 9m ln(l + e 

77 a =l *-JC 



/(«)(«) ^ 



(11) 



Due to the nonlinear nature of (g) and ([iT|), c e g(r) is usually obtained by solving the 
equations numerically. However, like the TBA equations, the NLIEs can be exactly evaluated 



at the ultraviolet point j28|, the result for the above equations being j25 | 



Ceff(0) = N 



3C 

e+i 



2 C ab la ««f» • 



(12) 



a,b=l 
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The effective central charge of the massive <3? a dj perturbation of WQ^ q is obtained by setting 
the parameter £ and the twists a = (a»i, a2, ■ ■ ■ , ajy) to 

Z = p/(q-p) , a = (2/p,2/p,...,2/p) . (13) 

The afflne Toda coupling constant [3 2 is related to £ via 1 = £/(£ + 1), and the above 
ensures fi 2 = p/q. The choice of a is motivated by (2^,|(]] for the A\ related models, and the 
prescription given in [25] which yields the central charge of the WQ minimal models rather 
than Cefj(r). As a first check we insert (|i~3l) into (12), simplify using 12 J2ab=i ^ab = Nh(h+1), 
and recover the expected UV effective central charge 

In accordance with (||) we modify the massive equations to interpolate from a model with 
£ = p/{q—p) to one with £' = (2p—q)/(p—q), that is £' = £—1. Motivated by |[],||], we associate 
two functions and to each node of the Dynkin diagram, introduce new kernels and 
twists such that the functions satisfy 

>(0) = 



JV 
6=1 



IV 



6=1 



d0'. 



C 'a6 la 6 



'a6l 



+Jd6' X ab( 



6>')3rmln(l + e 



(14) 



,(m) 



V 



(^) 



F 

6=1 



A" 



-2*E 

6=1 



/ d0' <t> ab (6-0')Sm ln(l + e-^^')) +/d0' x b(^-^) Qmln ( 1 + e/ * >( 
7c 7c 



and replace the formula for the effective central charge with 



Ceff(r-) 



R N 

%^2m a [ d9e e Qmln(l + e^ ?W) 

n a=l ^ C 



dee~ e Qmln(l + e 



(15) 



(16) 



As explained in 0, we fix the kernel functions by considering the equations in the limits 
in which r — > and r — > oo. In the far infrared, the massless equations coincide with the 
ultraviolet limit of the massive equations, with kernel </? o &(0) replaced by <t> a b{6)- Since this 
should describe a model with parameter £—1 we set cj> a b(0) to 

sinh(^fc) 



dk 



mt 



sinh(f (f-l)fc) cosh(p) 



(17) 



For very small r, the massless equations should instead coincide with the ultraviolet limit of 
the massive equations with parameter £. After some manipulations Jl|, the fourier transformed 
massive and massless equations can be directly compared. With the tilde denoting the fourier 
transformed functions, the equations will match provided Xab{k) and a' a satisfy 



N 



Vab{k) = (j) ab (k) + £ Xac(k)[l ~ Hk)]J 
c,d=l 

N 

"a = a' a + X: X«c(0)[l - 0(0)]^ Xdb(0) . 



(18) 
(19) 



c,d=l 
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Inserting the expressions for ipaj, and 4> a b into (|l8|) we find 

U H 'J-ml M w ) 2 °^ k) - £ ^ (fc> c " (t) "< fc> • (20) 

After multiplying both sides by Cf a (k) and summing over the index a, we can write the 
right hand side as a square: (x(k) C(k))'j b . Taking the square root and inverting the fourier 
transform yields 



= ± /^e^ ^^^rr^^^ ^w • (2i) 



dfe jke sinh(%fc) 1( 

2?r sinh(f (f-l)fc)cosh(f fc) afe 

The above only fixes Xab(&) up to a sign, but we find that choosing the negative sign results 
in an effective central charge consistent with a WQ minimal model, whereas the other choice 
does not yield a recognisable formula for c e ff (0). With the negative sign we find the new twists 
should be 

a' a = jzi aa ■ ( 22 ) 

To avoid the pole in (|2^) at £ = 1, and the poles in the kernels which cross the real axis as £ 
falls below one, we only consider models with 2p > q and therefore £ > 1. This is a sensible 
restriction since a flow of the form (||) with 2p < q would have an infrared CFT labelled by 
(2p — q,p), the first of which is negative. 

3. The massless NLIEs have ultraviolet and infrared values of c e ff(r) which exactly match 
those of the conjectured flow (||) provided we continue to use the massive prescription for £ 
and a @. We find 

, (23) 

The massless flows naturally fall into families indexed by an integer J = q—p : 

wg p / +J + $ adi ^wg p N - J * . (24) 

At J = 1 there is a unique family corresponding to the flows between the unitary minimal 
models. For these models we tested the massless NLIEs against the the TBA equations 
|1C|, typically finding very good agreement. As expected the effective central charge was 



consistently monotonic. For J > 1 there are f{J) different families, each of which interpolates 
between nonunitary models (here ip denotes the Euler-93 function). Solving the NLIEs for the 
nonunitary models, we found that c e g(r) increases away from its UV value, undergoes a 
number of oscillations and then settles down to the predicted IR value. The nonmonotonic 
behaviour of two families of flows is illustrated in Figures I andg 

4. To strengthen the validity of the conjectured equations we extract a number of further 
predictions, which are then compared with results from ultraviolet and infrared conformal 
perturbation theory. 

The UV groundstate energy of a CFT perturbed by a primary field $ of conformal dimen- 
sion Auv is predicted to behave as pi"|,|32| 

E (R) = RB(X) - ™igp , c pert (r) = c eS (0) + Y,C 3 (XR y ) j , (25) 
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Figure 111 A family based on WA3 with J=6. 
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Figure H: A family based on WD^ with J =2. 



where y = 2(1 — Atjv) and the coefficients Cj are proportional to the connected correlation 
functions of the perturbing field on the plane. The action (|6|) implies A and M must satisfy 

A = kM v (26) 

for a dimensionless constant k, further implying that c pe rt expands in powers of r y . On the 
other hand, the 2(£+l)i7r//i periodicity f\ of the nonlinear integral equations suggests that 
Ceff(^) expands as a series in r 2h ^ l+ ^\ which will agree with c per t(r) provided 2(1 — Atjv) = 
2/i/ (l+£). Substituting pj (q—p) for £, we find the NLIEs predict a value of Atjv which exactly 
matches the conformal dimension of <I> a( ij (||). 
The expansion 

00 

Ceff(r) = c eff (0) + B(r) + ]T c j W > ( 27 ) 

i=i 

differs from c pe rt( r ) by the bulk term B(r), but it may easily be extracted from the NLIEs. 
For this we need the leading asymptotics of the kernels as 9 — > — 00. The denominator of the 
inverse deformed Cartan matrix has either a term of the form sinh(7r/c) (A n ) or cosh(7rfc/2) 
(D n , Eq, £7, E%), both possibilities leading to a simple pole in the kernels 4> a b(9) and Xab{9) 
at k = —i. Moreover, the kernels (|l7| , ^l| ) have a pole at k = — i/t/(£— 1), which contributes a 
term proportional to e 11 ®/^ 1 ^. Therefore we have 

<M*) ~ ^ e e + p« e hd '^ + . . . , 9^ -00 , (28) 
X a& W~xil ) e g + cS ) e^/«- 1 ) + ... i 0--oo. (29) 

* While the periodicity is easily extracted from the associated Bethe ansatz equations J25|,[26), it is less trivial 
to see directly from the NLIEs. 
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The expansion coefficients can be calculated by taking the appropriate residues, but we will 
only need (j)^ and Xil > both of which can be written as 

AD _ ^(KWf) m AD _ sin2 (f) m {m) 

where 

!/(^-i) = l , !/(£>„) = 1/2 , = y273sin(^/12) 

^7) = 4 sin ( 7r / 18 ) sin (W9) , K^s) = - S ^ {n/m l, ■ (31) 
y3 2v3sm(7r/5) 

A small generalisation of the argument described in yields the massive and massless bulk 
terms: 

B It)- 3sin ^Mg) r2 fl (r) _ 3Kg) 2 

-"massive l,' J — n ■ /7r/i , c\\ ■ ' -"massless \> ) — . , n . 1 ■ l oz 7 

27rsm(f (1 + 0) sin(j) 2vrsm( I (l + £))) 

The bulk term corresponding to the massive perturbation of the unitary minimal models 
has been calculated by Fateev [33 in the context of the associated coset. By analytically 



continuing the coset parameter k to rational values (set k = £—h), we find Fateev's massive 
bulk term exactly coincides with ours for both unitary and nonunitary models. 

The massless bulk terms are new, but we can make at least one concrete check before 
turning to numerics: the models with p = h+1 and q = h+2 correspond to the coset Qh x 
Gn/G%, for which the massive and massless perturbations are known to coincide [|,34], as 
fortunately do our bulk terms at £ = h+1. 

The bulk terms have a simple pole whenever £+1 = mh for some integer m, which should 
cancel against one of the terms in the infinite sum so that the expansion ( p7|) continues to be 
regular 1 35 — 37 1 . In all cases the result, found by evaluating 

Jim ^f,- 2 -^], (33) 

is a logarithmic term: 

5 massiv e(r)| 5=m/t _i = — = r 2 lnr , 5 massless (r)| ?=m?t _i = r 2 lnr. (34) 

Now we are in position to compare (^5|) with (^7|) , apart from one remaining difficulty The 
perturbative coefficients Cj are usually hard to calculate, and instead it is easier to estimate, 
from the massless and massive NLIEs respectively, the expansion coefficients Cj and Cj. The 
perturbative coefficients Cj do not depend on A, and if we assume the mass and crossover 
scales M are equal we should find 

( i Ho- (35) 

We include a small sample of our numerics in table ^. Such good agreement provides an 
excellent check on both the massive and massless NLIEs, the associated bulk terms and the 
above assumption on M. 

5. Close to the infrared fixed point the model is described by the action of the infrared CFT 
plus an infinite number of contributions from irrelevant operators, resulting in a theory which 
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j 


c i 


C 3 


j 


C 3 


c j 





1.657142856 


1.657142856 





3.34117647 


3.34117646 


1 


-1.49526585 


1.49526587 


1 


-0.1962515 


0.1962514 


2 


0.01009599 


0.01009595 


2 


0.067572 


0.067573 


3 


-0.0015271 


0.0015278 


3 


0.003602 


-0.003603 


4 


0.0000603 


0.0000602 


4 


0.00005 


0.000008 



Table 2: Comparison of the massive and massless UV coefficients, found via the NLIEs. 



is unrenormalisable. However by considering the contribution of a finite number of fields it is 
still possible to make a comparison with results from either NLIEs such as ours or from TBA 



equations [38,3S]. We consider 

A = A m + g I ipd z x + t I TTd z x + ... , (36) 



where the (possibly missing) irrelevant field tp of dimension Air and TT of dimension 2 belong 
to the infrared CFT. The couplings are related to the crossover scale M as 

g = K g M 2 - 2AlR , t = K t M~ 2 . (37) 

The action implies c e g(r) has IR expansion 

oo oo 

Ccff(r) ~ c cfr (oo) + X>(M (2 ~ 2AmW + Y, f ^r)- 2j + ... . (38) 

3=1 3=1 

From the NLIE point of view Q corrections to c e ff (r) come from the 9 — ► — oo asymptotic of 
Xab{6) given by (|29|), the e h(> ' term generating a a series of the form r~ 2 ' l// ^ _1 ). Comparing 



with the CPT expansion (pq) leads to the prediction Air = 1 + h/(£—l), which can be 
identified with the conformal dimension of the primary field ^adj ' with weights O = Q^, O' = 

We can also extract a prediction for K t from the NLIEs. We start with the first two 
coefficients of the series generated by TT |39| : 



7r 3 C efr (oo) 2 7r 6 C e ff(oo) 3 

h = , h = , (39) 

and compare them with the coefficients of r~ 2 and r -4 found from the NLIEs. Adapting the 
TBA argument |i7]], we use the e e term in the expansion of Xab(6) to find 



Ceff(r) ~ Ce ff (oo) - ^c 2 qS (oo)^^ 2 + 2 C 3 ff (oo) (^xffl r" 4 + . . . . (40) 
The explicit form 



v« r- 2 = Sin2( ^ } r- 2 (41) 

Xlir 7rsin(^^-l))Ka) 1 j 
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indicates a pole whenever £— 1 = m'h for some integer m! . Evaluating as for the UV case we 
find ( pi] ) becomes 

(i) -2, 2(-ir , sin 2 (f) 
Xii r k=m'h+i = 2 TTr\^ T lnr • 42 

The infrared expansion coefficients have been obtained numerically for the models WA P,P+ + 
013 , p = 5, . . . , 10 in pp,[4p, and they show good agreement with our predictions, while for 
Wj\P'P+ anc j WD^^ ' n we find agreement with the theoretical results of [37]. Finally, the 



effective central charge of the </>i2,02i and 0i5 perturbations of the Virasoro minimal models 
is exactly half that of certain WA2 models: 

WA p ' q + cf> 12 <- WAfo + S^ 

WA™ + 4> 21 <- WAf p + <f> adi (43) 
WA^ + fa ~ WA^ + ^ adi . 

Only the (j) 2 \ and 0x5 perturbed models have a massless flow, and the IR expansion coeffi- 
cients found in Q match with those predicted above for the associated WA 2 model. The 
correspondence actually works for any value of £ since the 0i2/</>2i/0i5 NLIE is based on the 
tadpole diagram Ti, which is related by folding to A 2 . 

Finally, comparing (40, 41) to (38, ^9|) yields the promised prediction for nt- 



4sin 2 (f) 

Kt = 7r2sin(f(£-l)MS) • (44) 

6. We have shown that the function c c ff (r) for nonunitary models interpolates from a CFT 
in the ultraviolet to an infrared CFT that has smaller effective central charge, but in a 
nonmonotonic way. It is likely that there is a function, as yet unknown, which monotically 
interpolates between the nonunitary CFTs and satisfies a 'nonunitary c-theorem'. 

We'd like to make two further comments concerning the massive nonlinear integral equa- 



tions. First, level-rank (or KNS) duality [23,41] relates two nonunitary minimal models 



WA n n l 1 = WA p ~_^l 1 , n,pcoprime, (45) 

which both have a field <£ a dj with conformal dimension A = (pn—p—n 2 )/p, resulting in entirely 
equivalent perturbed theories. It is not at all obvious that the NLIEs based on A n _i and those 
based on A p _ n _i at the appropriate values of £ will produce the same value of c e ff (r), but 
provided the normalisation of the lightest mass is such that the coupling A is the same for 
both models, our numerical studies confirm this. By equating the massive bulk terms we find 

sin(vr/(p - n)) M [n _ x] = sin(7r/n) M [p _ ri _ 1] ; (46) 

the same result can also be deduced from the relation A = kM v given for the unitary minimal 
models in [33|. We can rule out a massless flow from WA^ P X via <I> a dj by considering the 



conjectured IR model W ' A 2 ™~y' n ', a sensible CFT if 2n—p > h. Since h = n this would require 
p < n, contradicting the assumption n < p. 

As mentioned above, the massive ^15 and 4>i 2 perturbations of the the Virasoro models 
WA X ' P (for odd p greater than 4) are integrable, and the groundstate energy can be found 
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using a single NLIE based on A^' 42], with appropriate tuning of the free parameters. Thus 
by level-rank duality the models WA^2^ P have up to two extra massive perturbations in 
addition to $ a dj whose groundstate energy can be described in terms of a nonlinear integral 
equation. 

Second, the massive NLIEs encode the finite-size effects of a further set of models obtained 
by sending £ = k—h to infinity. In this limit the perturbed coset x / + ^adj becomes 



the Q Gross Neveu model [22,43], with groundstate energy described in terms of an infinite 
number (k.N) of coupled TBA equations. The NLIEs offer a clear advantage over the TBA 
as the TV" equations can still be solved numerically at any value of r. The kernels 

S* e \ s *-^m c * (k) ) (47) 



form part of the prefactors of the associated S-matrices 44]. The perturbation is (almost) 



marginal and the ultraviolet expansion of the effective central charge no longer has a simple 



power series form (27). By studying the NLIEs (numerically and analytically) we hope to 
uncover the expected logarithmic corrections to c e fr(r). If this approach is successful there 
are a number of two-dimensional sigma models which have an interpretation as a perturbed 
conformal field theory [Hj. Since in all cases the groundstate energy is described by an infinite 
number of TBA equations, it would also be interesting to extend the current set of NLIEs to 
include such models. 

Finally we remark that even though the groundstate energy of the perturbed nonunitary 

(2) 

models is apparently real, recent results based on the A 2 models indicate the massive finite 



size spectrum may in general be complex f45| . It remains an open question to study the 
excited states of the massless models, perhaps via the massless NLIEs, to see if a similar 
result holds. 
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